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Abstract. The observations of the total electronic 
component and the positron fraction in cosmic rays 
by the FERMI, HESS, ATIC, and PAMELA in- 
struments are studied with analytical propagation 
models, both for a set of discrete sources and for 
a spatially smooth source distribution. The positron 
fraction over the entire energy range of 1 
100 GeV is shown to fit with the nested leaky 
box model. We derive the spectrum of electrons in 
cosmic rays arising from direct acceleration by the 
sources and discuss the narrow spectral feature in 
the spectrum. 
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I. Introduction 

Progressive improvements of detectors since the pi- 
oneering efforts in the 1960s have led to improved 
measurements of the spectrum of the electronic compo- 
nent and positi-on fraction by the FERMI, HESS, ATIC, 
and PAMELA instruments [1-5]. These measurements 
have stimulated great excitement, not only because of 
their importance to cosmic ray studies, but also because 
of the possible interpretation of some of the features 
from signals from the annihilation of dark matter in 
the Galaxy. We reproduce the observed total electronic 
spectrum ft{E) and the positron fraction R{E) in fig.[T] 
and fig. |2] respectively. 

In section II, we show that these observations help 
us in distinguishing between the leaky box ||6l and the 
nested leaky box models Q. It is the difference in the 
kinematics of the reactions producing secondary 'light' 
nuclei and positrons that helps in distinguishing amongst 
these models. We show that the positron fraction ob- 
served by PAMELA may be fit well with the nested 
leaky box model. We begin section III by deriving the 
spectrum of electrons resulting from direct acceleration 
in the cosmic ray sources by subtracting the secondary 
electrons and positrons from the spectrum of the total 
electronic component. After this, the rest of the section is 
devoted to a discussion of this primary electron spectrum 
in terms of a set of discrete sources of cosmic rays 
sprinkled throughout the Galaxy. The analysis indicates 
that we need the closest source to be no more than ^ 
200 pc from the solar system, which in turn implies 
that there are about 5000 sources accelerating cosmic 
rays in the Galaxy. The contrast between the spectra 
expected from a discrete set of sources and from sources 
spatially distributed smoothly throughout the Galaxy 
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Fig. 1 . The measurements of the total electron spectrum, a smooth 
fit (solid line), and the total secondary electron and positron spectrum 
calculated from the nested leaky box model (thick solid and dotted 
lines) for Ec = 600 GeV and Ec = 1200 GeV. 
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Fig. 2. The theoretically calculated positron fraction in models A 
(similar to that of Moskalenko and Strong |T0|), B, and C are compared 
with the observations. All calculations are normalized at ~ 10 GeV. 

plays an important role in analyzing signatures from the 
annihilation of dark matter in the Galaxy. This aspect 
is discussed in section IV. It would be appropriate to 
note here that the theoretical considerations of cosmic 
ray transport with exponential pathlength distributions 
|6||8| and its analysis in terms of a set of discrete 
sources sprinkled over the galaxy HH) started three or 
four decades ago. The recent improved measurements 
constitute a testing ground for these early ideas on 
cosmic ray physics. 

II. Positron Fraction in Cosmic Rays and the 
Leaky Box Models 

The observations of decreasing ratios of secondaries 
to primaries like B/C and (VH-TiH-Sc)/Fe with increasing 
energy per nucleon in cosmic rays could be incorporated 
into the original leaky box model (A and B) |6| by 
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Fig. 3. The observed B/C secondary to primary ratio is plotted (points 
from a compilation in |11J) along with the power law extrapolation 
at high energies (dot-dashed line. Model A), a constant extrapolation 
(solid line, Model B), and a two-component fit (dotted lines. Model 
C). 



letting the mean residence time r depend on energy 
in such a way as to reproduce the observations [11]. 
A viable alternative to this idea was the nested leaky 
box model (C) Q in which the energy dependant part 
of the ratio (proportional to Ts{E)) was attributed to 
the spallation in the matter constrained in a cocoon- 
like region surrounding the sources and the resident 
time in the Galaxy tq was assumed independent of 
energy. With observations of the ratios limited to a finite 
range of energies, it became customary to extrapolate 
the residence times as a simple power law, decreasing 
as ~ tE^^. Even though the observations to date allow 
A beyond some tens of GeV/n, this alternative had 
not been discussed earlier. These three possibilities and 
their corresponding fit to the secondary to primary ratios 
are shown in fig. |3]for the B/C ratio. The (V-i-Ti-i-Sc)/Fe 
ratio can be fit similarly. The three models presented 
here are characterized by mean residence times as noted 
below: 

A: ta{E) - tqE--^, E>2 GeV/n (1) 

TBiE) ^ta{E), E<lQGeV/n 
^TG^ X, E>10 GeV/n 

(2) 

Ts{E) ^tb{E)-tg, E<\()GeVln 

^TG^X, E>10 GeV/n ^ ' 

where X is a constant. 

Consider now the injection of cosmic rays by sources 
with the spectral form 

Q{E) - QoE^'^cm^^s'^sr~^GeV/n'^. (4) 

When spallation and energy loss during transport could 
be neglected, the three models led to the following 
equilibrium spectra of cosmic rays in the Galaxy: 



A:fA ^ Q^T,E-if'+^) 

B: fB - QQTQE-''fi+^\ E<IQ GeV 

- QqtgE-P, E>\Q GeV 
C: fc - Q^tgE-P. 

Note that fc does not depend upon Ts{E) Q. 



(5) 



Unless /3 itself is a function of energy so that eq. 5B 
reproduces the observed spectra of the nuclear compo- 
nent that is a simple power law, model B is not viable. 
Model A requires /? ^ 2.22 and predicts increasing 
anisotropy with energy of cosmic rays, and model C 
requires (3 ~ 2.65 and generates constant anisotropy at 
all energies. 

The observations of the positron fraction provides a 
clear way of choosing amongst these three models. The 
reason for this is related to the fact that in the spallation 
process that generates the secondary cosmic ray nuclei 
like B in collisions of the primary C nuclei, the daughter 
nuclei will emerge with the same energy per nucleon 
as the parents. On the other hand, the production of 
positrons proceeds through the production of mesons 
(mainly pions) in the collision of the primary nuclei 



which follow the decay chain tt^ 



followed by /i 



process, the positrons and the secondary electrons carry, 
on the average, only a small fraction ~ 0.05 of the 
energy/nucleon of the primary. As a consequence of this, 
the production spectrum of positrons (and secondary 
electrons) in all the three models are nearly identical. 
The fact that low energy cosmic rays spend more time 
in the cocoon surrounding the sources plays no role 
as the e+ and e~ secondaries are generated only by 
the high energy part of the nucleon spectrum. Now the 
equilibrium spectrum of secondary e+ and predicted 
by the three models simply follow the product of the 
production rate, Qe±{E) and the residence time t{E) 
at low energies. At very high energies, all the models 
predict steeper spectra ~ Qe±{E)E^^: 

A: fn±{E) = Qn±iE)TAiE) 

^ ToE-^^+^\ E^E, 

B: fn±{E) = Qn±{E)TB{E) 

- E<lQGeV 

- E-f^, 10 GeV <E^Ea 

- E-'^+'^, E:$>Ec 

C:fn±{E) = Qn±{E)TG{E) 

- E-f\ E<t:E^ 

- E-^^+\ E:^Ec (6) 

where Ec ~ l/brG and b is defined in eq. 10. In fig. |2] 
we show the positron fraction 



R{E) = U+{E)/ft{E) 



(7) 



with the positron spectrum /„+(£') as given in eq. 6 for 
the three models, and the spectrum of the total electronic 
component is taken to be the smooth fit to the data shown 
in fig - [n Note that only model C provides a good fit over 
the entire range of the PAMELA observations. Thus, 
the observed positron fraction resolves the degeneracy 
amongst models developed to understand the nuclear 
secondaries and suggests that the nested leaky box 



PROCEEDINGS OF THE 31"* ICRC, LODZ 2009 



3 



model may provide a closer approximation to cosmic 
ray transport in the Galaxy lfT2l . 

III. Spectrum of Electrons Generated by 
Cosmic Ray Sources 

The nucleonic component in the primary cosmic rays 
consists mainly of protons with a small fraction of 
neutrons which come in bound to He and other nuclei. 
As a consequence, the production of positrons is favored 
over that of electrons in cosmic ray secondaries: 

77 = Q„_(£;)/Q„+(S)«0.5-0.8. (8) 

The lower value of is favored by the theoretical esti- 
mates based on the data on nuclear interactions obtained 
at accelerators [131, and the larger value is favored by 
the observations of the iJ,~ ratio generated by cosmic 
rays in the atmosphere [\A\. The loss of energy by 
energetic cosmic ray electrons due to inverse Compton 
scattering and synchrotron emission in the magnetic 
fields prevalent in the propagation volume is given by 

dE/dt = -bE^; (9) 

b = (3.22X10-V— 

V \eV ■ cm / 

+7.9 X 10-'Hlaauss)GeV-H4yr-\(lO) 

Noting that the scattering of starlight by the very high 
energy electrons is described by the Klein-Nishina for- 
mula, the effect of starlight may be neglected, and only 
the microwave background contributes to ujph in eq. 10. 
Taking H ^ 5 fiGauss and u)ph ~ 0.25eV cm^^, we 
get b « 1.56 GeV~^Myr~^ . Adopting the nested leaky 
box model, the spectrum of the secondary positrons and 
electrons is given by ID 

U+{E) = fn-{E)/r, (11) 

pl/bE 

= / QoE-^{l-hEtf-^e-'/^'^dt. 
Jo 

In fig. [U we show the total secondary component 
fs{E) — fn+{E){l + 77), with the normalization for 
Qo determined from a smooth fit to the PAMELA data 
and TG ~ IMyr, i.e. E^ = I /bra « 600 GeV. We also 
show the secondary component for Ec ~ 1200 GeV . 
Subtracting this (for Ec — 600Gew) from the total spec- 
trum of the electronic component, we get the spectrum 
due to the sources, 

gc{E)^ ft{E)- U+{E){l+T^), (12) 

which is displayed in figs. |4] and |5] 

In an attempt to understand this spectrum, we assume 
diffusive transport of cosmic rays and write the differ- 
ential equation describing the transport as 

dN/dt- nV^N + N/t = Q. (13) 

Note that here we have simplified the form of the 
Green's function by introducing an escape term t/r 
in place of the boundary condition demanding that the 
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Fig. 4. The primary electron spectra due to a single source at various 
distances from the source with Ex = 5 TeV compared to the primary 
electron spectrum. 
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Fig. 5. The primary electron spectra resulting from many cosmic 
ray sources at various distances to the first source with Ex = 5 TeV 
compared with the primary electron spectrum. Their difference (thick 
line) is shown. 

cosmic ray density vanish at the planar surfaces of 
the thick cosmic ray disk of the Galaxy. The Green's 
function is given by 

G(r, t) = {'iT:nty^''^exp(^ ~ ~ 7) " 

This along with the subsidiary eq. 9 is adequate to 
describe the transport of the electronic component of 
cosmic rays. Note that for a smooth and uniform distri- 
bution of sources Q{r,t) — S{t), the Green's function 
in eq. 14 integrates to the simple leaky box model. 

Suppose we have a source situated at a distance 
which is continuously emitting cosmic ray electrons with 
a spectrum 

Q ^ QqE~^, E < Ex 

^ 0, E>Ex, ^'"^ 

the observed spectrum is given by 

QoE-^{l - bEtf-^Gin, t)dt. 

(16) 

We display in fig. |4] this spectrum for various values 
of r,. Note that fov T E^ Ki b - 10 TeV, and 

ri K, 0.2 — 0.5 kpc, one can obtain a reasonable fit to 
the high energy part of the primary electron spectrum 
g(E). 
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Fig. 6. The spectrum of electrons accelerated in shocks to fit the 
electron excess for a cutoff energy between 5 — 10 TeV (thin lines). 
Also plotted is the spectrum for a (5-function input from continuous 
distribution of sources with Ea = 1200 GeV. Note the peak at ~ 
300 GeV ^ Ec/2 (thick dot-dashed hne). 
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Fig. 7. Shown here is the spectrum from a (5-function input from a 
single source (solid lines) and the spectra from the sum over a discrete 
distribution of <5-function sources (dotted lines) for various distances. 



However, it is appropriate to add the contributions of 
all the sources situated at different distances to obtain 
their net contribution foiE) to ge{E). Assuming that 
the i*'* source is at an average distance of Viri, we 
sum their contributions and compare with the primary 
electron spectrum ge{E) in fig.|5] We see that irrespec- 
tive of the choice of parameters, the high energy part of 
ge{E) is not reproduced by the theoretical expectation 
foiE). This difference becomes particularly large for 
ri much larger than ~ 200 pc. Thus we may expect that 
the distance to the nearest source and the typical spacing 
between the sources is ^ 200 pc. The difference between 
the contribution of the discrete sources foiE) and the 
fit to ge{E) is also shown in fig. |5] and the difference 
ne{E) with the data points is shown in fig. |6] 

IV. Narrow Spectral Features in the Primary 
Electron Spectrum 

The narrow spectral feature displayed in figs. |5] and 
|6] has been ascribed in toto or in part to products of 
dark matter annihilation (see ifTSl for references). In this 
section, we discuss two input spectra, one a (5-function 
in energy and the other a flat spectrum ^ E^"^, such as 
that expected for acceleration at planar shocks of high 
Mach number. 



A. 5-function Input Spectrum 

Here, the injection spectrum Q{t = 0,E), the spec- 
trum after time t, fsi{E) (without leakage), and fsc{E) 
for a spatially smooth source and leakage lifetime tq 
are given by 



hEt 

En - 



'E, 
E 



(17) 
(18) 



V 



(19) 



bE^ KbEaETcJ' 
The spectrum fsc{E) peaks at £' = Ec/1 = l/bTc for 
Ea > Ec/2, and there is no sharp peak at Ea (see fig. 
|6]l. A (5-function input from discrete sources leads to an 
observed spectrum 

^ ,^ , E^af bEaE \3/2 



.4TTK{Ea'E) 

brfEEn 



xexp- 



V4k 



Ea-E 



i{Ea-E) ' bEaEra)' 
For large tq, this spectrum displays a peak at 

GnEa 



E.^ 



'peak 



(20) 



(21) 



6k + br'^Ea ■ 

For small r^, the peak will be sharp near Ea, but with 
increasing r^, the peak will shift to lower energies and 
will become broader. We display in fig. |7] some examples 
of the spectra generated by (5-function inputs at various 
Ti and show the sum, fo{E), over the sources at various 
distances in fig. |5] 

B. Shock Acceleration 

Planar shocks of high Mach number yield 

Qshock{E) ^ QoE ^, E < Ex 

resulting in an equilibrium spectrum 



(22) 



f2{E) 



E^ 



1 



(23) 



which can reproduce ne{E) as shown in fig. |6] 

We refer the reader to ref. [15] for a comprehensive 
overview of current efforts to explain these observations 
and to ref. [16] for more details regarding this paper. 
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